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Dynamical heterogeneity is considered to be key in the understanding of the dynamical behavior 
of slowly relaxing disordered materials. In this work, we test one hypothesis for its origin, namely 
that it emerges from soft (Goldstone) modes associated with a broken continuous symmetry under 
time reparametrizations. We do this by constructing coarse grained observables and decomposing 
the fluctuations of these observables into transverse components, which are associated with the pos- 
tulated time-fluctuation soft modes, and a longitudinal component, which represents the rest of the 
fluctuations. Our test is performed on data obtained in simulations of four models of structural 
glasses. As the hypothesis predicts, we find that the time reparametrization fluctuations become 
increasingly dominant as temperature is lowered and timescales are increased. Especially, the ra- 
tio between the strengths of the transverse fluctuations and the longitudinal fluctuations grows 
as a function of the dynamical susceptibility, X4: which represents the strength of the dynamical 
heterogeneity; and the correlation volumes for the transverse fluctuations are approximately propor- 
tional to those for the dynamical heterogeneity, while the correlation volumes for the longitudinal 
fluctuations remain small and approximately constant. 



I. INTRODUCTION 

The rapidly increasing relaxation timescales, the presence of non-exponential relaxation, as well as the violation of 
Stokes-Einstein relations between viscosity and diffusivity are some features observed close to the glass transition in 
glassy systems [TH3]- The appearance of these features suggests that relaxation dynamics is heterogeneous, i.e. that 
it is faster in some regions and slower in others [2H11- Direct microscopic evidence for this behavior has been found 
both in simulations [H |5] and in experiments [H [6Hl6] . The understanding of dynamical heterogeneity is believed to 
be crucial to explain anomalous behavior of materials near the glass transition, and even possibly to explain the very 
presence of the glass transition itself (2| . Despite many efforts trying to address the origin of dynamical heterogeneity, 
this question still remains open [H [T7lf20] . 

In recent years, different tools have been used to probe dynamical heterogeneity. One such tool is the dynamical 
susceptibility, X4, which depends both on the strength and the spatial extent of the fluctuations. The peak value of 
X4 has been observed to grow while approaching the glass transition. However, this tool by itself cannot tell us much 
about the origin of dynamical heterogeneity and it may be desirable to supplement it with other ways of probing the 
fluctuations in the system. 

In this work, we test the predictions of a theoretical framework that aims to describe the slow part of the fluctuations 
in the relaxation dynamics |20H24) . This framework is based on the hypothesis that in glassy systems the long time 
dynamics is invariant under reparametrizations of the time variable |20[ I21|, l25H28j , but this invariance is broken, 
giving rise to Goldstone modes that manifest themselves with the emerge of heterogeneous dynamics [5DH23J • The 
Goldstone modes correspond to fluctuations in the time reparametrization t — >■ 4>rit), i.e. 

Cp{t,t')^c[Mt),Mt% (1) 

where C{t,t') is a global two-time correlation function. Some indirect evidence in favor of the presence of this kind 
of fluctuation in structural glasses has been presented in ;29n32j. In this work, we present results for a more direct 
test, based on decomposing fluctuations into a transverse part satisfying Eq. (1) and a longitudinal part containing 
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all other fluctuations. This procedure allows one to separately quantify the strength and spatial correlations of 
both kinds of fluctuations, as a function of temperature and timescales, for a variety of glass-forming models both 
below and above the mode coupling critical temperature T^. The same kind of test can also be applied to particle 
tracking experimental data from colloidal and granular systems [^I5H15|. thus allowing to investigate a possible unified 
explanation of dynamical heterogeneity in diverse systems. A summary of the results of an early version of our analysis 
was published in Ref. [33| . 

This manuscript is organized in the following way. In Sec. [IT] we discuss the hypothesis and define the quantities we 
use to test it. In Sec. |III| we present the details of the numerical simulations used to test the hypothesis and present 



the results of our analysis. Finally, we discussed our conclusions in Sec. IV 



II. TIME REPARAMETRIZATION INVARIANCE AND FLUCTUATIONS 



We start by discussing in general terms the expected effects of time reparametrization invariance on the fluctuations. 
Let's consider the global two-time correlation function 



1 ^ 

^(^'*') = ]^ E^o^i'?'- m) - rj{t')]). (2) 

where N is the number of particles, rj{t) is the position of particle j at time t and g is the wave- vector that corresponds 
to the main peak of the static structure factor. In equilibrium, and more generally in time translationally invariant 
(TTI) systems, this correlation function depends only on the difference between the two times t — t' . In the case of 
aging, the system is no longer TTI and C{t,t') depends nontrivially on both times t and t'. 
To analyze the fluctuations in the dynamics we define the local correlation as 



^rM^J^^ E C08{q-[r,{t)-r,{t')]}, (3) 

where the average over particles in Eq. ^ is now restricted to a region Bp around point r, which contains N{Bp) 
particles at time t' < t. Cp(t,t') probes the mobility of the particles in a region near r between times i' and t: it 
is close to zero for regions where the particle configuration has changed significantly, and much closer to unity for 
regions where it has changed little or not at all. 

In this work we will consider a slightly more restrictive version of Eq. ([T]), namely |33fl35| 

Cr-{t,t') ^ g[Mt) - Mt% (4) 

which corresponds to the case when the global correlation has the form [33H35] 

C{t,t')^gm)~^{t')]. (5) 

In principle, the functions g{x) in Eq. Q and Eq. ^ could be different. For example, it has been claimed [2j that non- 
exponential global relaxation could be the result of combining local exponential relaxations with different relaxation 
times, in which case we would have ^global = ^exp(— |a;|^) in Eq. ^ and giocai = A'exp(— jxl) in Eq. Q 'M'. In 
our case, for simplicity, we impose ggiobai(a^) = gioca.\{x) = g{x). Both this restriction and the restriction imposed in 
Eq. ([5]) could in principle make the results appear to be slightly less consistent with the hypothesis than they would 
be otherwise. The forms of the functions g{x) and can be obtained by fitting the global correlation C{t,t'), as 
we do in the next section. We will assume, for the moment, that g{x) and 4>{t) are known. 

As mentioned before, we refer to the fluctuations described by Eq. Q as transverse fluctuations and the rest of 
the fluctuations as longitudinal fluctuations |37j . To visualize these ideas in more detail, we define 

W=5"MCr(ta,ifc)], (6) 

where ^ab.r, with a, 6 € {1, 2, 3} fSS*, is composed of transverse fluctuations, 

^ab,r^ = Mta)- Mtb), (7) 

and longitudinal fluctuations, ^ab.r^, i-e. 

g~^[Cpita,th)] = $afa,r- = '^>ab,r^ + $afc,/ = M^a) " Mh) + $af,,/- (8) 
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This means that in the absence of longitudinal fluctuations, Eq. Q would be exact. In order to quantify both kinds 
of fluctuations, we define 

<yp = ^ {g~^[Cr{tiM)]+ 9~\Cp{h,h)]- g-^c^iM)]} 



V3 



^ ($12,/ + $23,/ - $13,/) , (9) 



= ^ [($12,/ + $12,/) - ($23,/ + $23,/)] (10) 



and 



712,,. = -^{5-MC,-(*i,t2)] + <?-MaK<2,i3)]+25-MCK*i,i3)]} 

= [($12,/ + $12,/) + ($23,/ + $23,/) + 2($13,/ + $13,/)] , (H) 

with t\ < ±2 < t^. As shown above, ap contains only longitudinal fluctuations, but tti^p and 7r2,f contain both 
transverse and longitudinal components. If Eq. Q was an exact identity, the local two-time function would verify 
the following relation 

9-'[Cp{h,t2)] + g-'[Cp{t2,h)]-g-'[Cpih,t3)] = 0. (12) 

Then, in the case where no longitudinal fluctuations are present, the vector {ap, tti ^, 7r2,f) would be restricted to be 
fluctuating in the plane ap = 0. We expect that as the temperature becomes lower, the timescales become longer, and 
the system becomes more glassy, transverse fluctuations should become progressively more dominant, according to the 
hypothesis (Eq. ([T])) and therefore the probability distribution p{ap,iTi^p,iT2,p) should become anisotropic, extending 
mostly along the ap = plane and not away from it. 

An analogous set of quantities can be defined starting from the global correlation C{t, t'). We can write the analogs 
of Eqs. ^ and ([8|, i.e. 

$a/ ='^(ia)-</'(ib), (13) 
g-^[C{ta,h)] = $ab = $a/ + $a/ = ^ta) ~~ (j^ih) + $a/. (14) 

However, in this case, the meanings of the symbols are different, ^ab^ and ^ab^ do not fiuctuate, because they are 
Dmputed from the global correlation in the thermodynamic limit, whic h is self averaging. If Eq. ^ was exact, then 



^ab = 0. However, in practice Eq. ([5|) is only approximate: in Sec. HI we fit the Ihs by an expression with the 
form given in the rhs. Thus ^ab^ represents the part of ^ab that can be represented as a difference (pita) — 4>{tb) 
according to the fits, and <I>ab^ represents the part of ^ab that the fit does not capture, or in other words a fitting 
error. Additionally, the global cr, tti and 1:2 can be defined, by analogy to Eqs. (|9|, ([To| and (11), in the following 
way: 

CT EE ($12 + $23 - $13) - ^ ($12^ + $23^ - $13^) , (15) 

TTi = ($12 - $23) , (16) 

and 
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TTs = 4^ ($12 + $23 + 2$13) , (17) 

with ti < t2 < t^. If Eq. ([5| was an exact identity, i.e. if the fit of C{t,t') by g[4>{t) — 4>{i')] was perfect, then ct = 
and the global two-time function would verify the relation 



9-^[C{hM)]+ 9-^[C{t2M)]- 9-'[C{tiM)] = 0. 



(18) 



A mor e ex tensi ve a nalysis of the fluctuations can be performed by separating longitudinal and transverse components 
in Eqs. (10) and (11). First, we recall that cr is a purely longitudinal quantity (see Eq. (joj)), therefore ap^ — and 



$13, f?-^ — $12, r"'" + $23, f"^- Now, wc make the following two assumptions. One is that the transverse and longitudinal 
fluctuations are not correlated to each other. The other is that all slow fluctuations are transverse, or in other 
words, that longitudinal fluctuations are short range correlated in time, or at least they are correlated over times that 
are shorter than the shortest time interval between the configurations that are being considered. This leads to the 
conditions 



(5$a(,,F^<5$cd,r-''^) =0 V a, 6, c, d, f, 



(19) 



and 



('5$a6,F^(S$c<i,r'^) = ioT a ^ c ov b ^ d, V 



(20) 



Here 5x = x — (x), and (■ • • ) denotes an average over thermal fluctuations, which in our case is approximated by 
an average over independent molecular dynamics (MD) runs. By using Eq. 0, it can be shown that the transverse 
components of Eqs. (10) and (111 are given by 



7ri,r 



;^($12,/-$23,r-'^) 



(21) 



and 



V6 



($12,r-'^ + $23,r^)- 



(22) 



Regarding the longitudinal components, by using Eq. (20) it can be shown that 



{(ScTf'^r) = ^ [((^$12,/)') + (('5<i>23,/)'> + ((^$13,/)')] 



(23) 



{(Snuf'r) = \ [(('5*12/)^) + ((^$23./)^)] 



(24) 



and 



((<5^2/)') = g [((<5$12,/)') + ((5$23/)') +4(((5$i3,/)2)] 

From the last three equations it can be easily seen that 

{{S^i/f + {S^2/f - 2(<5a/)2) = 0. 



(25) 



(26) 



By using Eq. ( 19 ) we find that, all together, we can compute the variance of the transverse fluctuations by combining 
the fluctuations of Eqs. (|9|, ( 10 ) and ( 11 ) in the following way 



(27) 
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Further, we can estimate correlation volumes by using the formula 



(28) 



where X4,aj5 = ^ '^r is a local coarse grained variable, a = J j^-ap, V is the volume of the system and Vcg is 

the coarse graining volume. Therefore, by using these equations together with Eq. (27) we can estimate the correlation 
volume of transverse and longitudinal fluctuations, respectively given by 



yT _ X4,Tri_i,T 



(29) 



and 



v{m?) 

Vcg {{bo?Y) Vcg ((<5a^.)2) ■ 



(30) 



If the time reparametrization hypothesis is correct, we expect that the variance as well as the correlation volume 
of the transverse fluctuations will grow together with those corresponding to the dynamical heterogeneities. 



III. TESTING THE HYPOTHESIS 



A. Systems 



We performed classical Molecular Dynamics simulations in systems of N particles that were initially equilibrated at 
high temperature Ti ^ Tc (where Tc is the Mode Coupling critical temperature [38,J, then instantaneously quenched to 
a final temperature T and allowed to evolve for times several orders of magnitude longer than their typical vibrational 
times We generated eight datasets by simulating four atomistic glass-forming models [52], which include 80:20 

mixtures of particles interacting via either Lennard- Jones (LJ) potentials 39J or via purely repulsive Weeks-Chandler- 
Andersen (WCA) potentials, as well as short (10- monomer) polymer systems interacting via either LJ potentials or 
via WCA potentials. The details of the simulations are summarized in Table |l] 









Systems 






Label 


N 


m 


Potential 




Runs 


A 


8000 


10 


LJ 


^ 0.9 


100 


B 


8000 


10 


WCA 


^ 0.9 


800 


C 


8000 


1 


LJ 


« 0.9 


250 


D 


1000 


1 


WCA 


« 0.9 


5000 


E 


1000 


1 


WCA 


« 1.1 


9000 


F 


1000 


1 


WCA - Eq 


« 1.1 


9000 


G 


1000 


1 


WCA 


« 1.5 


4999 


H 


1000 


1 


WCA - Eq 


^ 1.5 


4999 



TABLE L Details of the numerical simulations analyzed in this work. We considered systems of N particles, with 

TO particles per molecule, interacting via either Lennard- Jones (LJ) potentials or via purely repulsive 
Weeks-Chandler-Andersen (WCA) potentials, at final temperature T. Each temperature is described by its ratio 
with respect to the empirically determined Mode Coupling critical temperature Tc [38 for the same system. The 
last column lists the number of independent runs. Datasets F and H correspond to systems in equilibrium, all others 

to systems in the aging regime. 
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B. Results 



We begin by identifying the functions g{x) and (Eq. ([5])) that best describe the global correlations C{t,t') 
computed from our datasets by using Eq. We find that, for all datasets, the two-time correlation C{t,t') can be 
well fitted by using the form g{x) = qEA^xp [—{x/Oq)^]. However, as shown in the inset of Fig. (fll), the relaxation 
for different systems presents different behaviors, which leads to different forms of (f>(t). The best fits of C{t,t') that 
we obtained correspond to the following forms: for aging polymer systems (f)(t) = \n°'(t/to), for aging particle systems 
(j){t) = {t/to)", and for equilibrium (j>{t) — t/tQ. We can verify our proposed Eq. ^ for the different cases by using 
the form of the obtained functions g{x) and (j){t). For equilibrium we trivially recovered, as expected, the case of TTI. 
In the aging cases we can be verify that Eq. ^ can be rewritten in the form C{t,t') — C [h(t) /h{t')] [33], which is 
found in many aging systems |34j . Once the fitting procedure is performed for a given dataset, we can use the known 
values of the parameters Qea, P and to compute 

<^ab = 9-^[C{ta,h)] = eo{-\n[q~E\C{ta,h)]Y'^ (31) 

and, by using Eqs. (1517), to compute cr, tti, and 7:2- Fig. ([T]) shows the results of plotting the global values of a 
against the global values of ^tt^ + 7r| for all times, ti < t2 < ^3, and all systems. As discussed before, since the fits 
are not perfect, the results do not fall exactly on the line cr = 0, but the collapse and the fits are good enough to 
allow us to test the hypothesis. 
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FIG. 1: Global (cr, 7r|) pairs for all datasets (as labeled in the key) and all possible times ti <t2 < t^. Inset: 

Relaxation time r as a function of the waiting time tw for all the models and temperatures considered, as identified 
in the key to the main panel. r(tu,) is defined by the condition C'(tw + T{tw),tw) = 1/e. 



As mentioned in the previous section, if Eq. ([T]) is satisfied, we expect the local quantities to satisfy dp = 0. In 
Fig. ([2]) we show the results of plotting for all systems the 2D projection of the joint probability density of the coarse- 
grained local correlations, pdAcr^?], IAtt^I), with |Acrf?| = |crf— cr| and IAtt^I = ^J (j^i,? — tti)^ + (7r2,r — t^2)'^- The values 
of the global quantities are subtracted from the local quantities to avoid trivial effects due to differences in the global 
values. By doing this we are able to better compare the contours independently of the choice of C(ti, 12) and C(i2, ^3)- 
We do, however, keep the value of C(ti, t^) approximately the same for all systems, in this case C(ii, ^3) ~ 0.23. The 
three contours shown for each dataset enclose respectively 25%, 50% and 75% of the total probability. We coarse grain 
over moderately large regions, containing on average 125 particles, in order to detect collective fluctuations and, since 
time reparametrization symmetry is a long time asymptotic effect, we choose the times as late as possible. As the 
time reparametrization hypothesis predicts, for datasets with T/Tc ~ 0.9 (A-D), the purely longitudinal fluctuations 
Acrp are clearly smaller than the Att^? fluctuations, which contain both transverse and longitudinal contributions. This 
behavior is more noticeable for the 25% contour, which encloses the most likely fluctuations, than for the 50% and 
75% contours. For moderately higher temperature, T/Tc ~ 1.1 (datasets E and F), the anisotropy |Acrp| < |A7r^?| 
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FIG. 2: 2D contours of constant joint probability density p(\/S.ap\, computed using coarse graining boxes 

containing 125 particles on average. By subtracting the global quantities to the local quantities we avoid trivial 
effects due to differences in the global values. Each panel from A to H contains results from the corresponding 
dataset, for C(ii,i3) « 0.23, with the times chosen as late as possible within each dataset. Each set of three 
concentric contours is chosen so that they enclose 25%, 50% and 75% of the total probability. 



is still present, but less pronounced. In the case of the highest temperature, T/Tc ~ 1.5 (datasets G and H), the 
anisotropy is either very slight, or absent. In the case of the systems that are equilibrated, F and H, we find that the 
shapes of their contours are similar but slightly more anisotropic then the ones obtained for the same temperature in 
the aging regime, E and G, respectively. The effect of temperature in the anisotropy of the contours can be observed 
in more detail in Fig. ([3^) where the 25% contour of the probability density for the systems of particles with WCA 
interactions is shown for three temperatures (datasets D, E, G). This can be directly connected to the fact that, as 
the temperature is increased, the separation of timescales is less pronounced, the finite time corrections to the time 
reparametrization symmetry become larger, and the effects of local time variable fiuctuations become weaker. The 
same trends can be observed in Figs, ([sjs) and ^\p) for the same datasets as in Fig. ([3^) but for different values of 
the global correlation C{ti,t^). 

We can further analyze the effects of choosing different conditions from the ones chosen in Fig. ([2]), for instance, by 
comparing the results shown in Fig. ^ with results obtained for smaller course graining regions or for shorter times 
in the aging regime. Exactly this kind of comparison is shown in Fig. Q, where the 25% probability contours for 
dataset B are shown for three conditions. The contour labeled B is the one shown already in Fig. ([2]). The contour 
labeled B' corresponds to the same time, but with coarse graining regions containing on average 23 particles instead 
of 125. This leads to less averaging and stronger fluctuations, but also, since fluctuations correlated over shorter 
distances are no longer preferentially suppressed, the shape of the contour is no longer dominated by collective modes, 
and thus contour B' extends more in the direction of \^^<Jp\ than contour B. The contour labeled B" corresponds to 
the same coarse graining size of contour B, but with much shorter times. This leads to stronger finite time effects, 
analogous to the ones found at slightly higher temperatures, and as expected the contour is less anisotropic , and 
indeed, it resembles the contours corresponding to T/Tc ~ 1.1. 

We now move to a more quantitative analysis of the strength and spatial correlations of the fluctuations by making 
use of the results derived in Sec. |llj In the left panel of Fig. ^ we show the ratio between the variances of local 
transverse fluctuations and longitudinal fluctuations, (^{Stti^p^)'^ + [Sir^^f'^)'^) / {{Sapf^ ("variance anisotropy ratio") 
(see Eq. (27l), as a function of the initial time ti. Similarly, in the right panel of Fig. ^ we plot the ratio between the 
correlation volumes of transverse and longitudinal fluctuations, V"^ /V^ ("correlation volume anisotropy ratio") (see 
Eqs. (29l and (30)), also as a function of ti. We find that for aging systems both ratios grow as ti increases, as one 
could expect from the fact that at later times the reparametrizations symmetry breaking terms in the action should 
become progressively weaker |21| . For the equilibrium datasets, the dynamics is time translation invariant (TTI), and 
we observe, as expected, that both anisotropics are independent of ti. In Fig. ([6]), we show the same two ratios as 
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FIG. 3: Comparison of the contours enclosing 25% of the probabihty for the aging WCA systems at three different 
temperatures (D, E and G). The three different panels correspond to different values of C{ti,t^). 



functions of the strength of the dynamical heterogeneities, measured by X4 [C(ii, is)] = N (C(ii, ia)^) — {C{ti,t^)Y 

We observe that both anisotropy ratios grow when xa [C{ti^ t^)] increases, i.e. as the dynamical heterogeneity becomes 
more pronounced. Although the same qualitative behavior is observed for all datasets, the curves are different for 
different systems and temperatures. 

The results presented in Fig. ^ correspond to the value C(ii,i3) ~ 0.23, but similar results can be obtain for 
different values of (7(^1,^3), as one can expect from Fig. Results for three different values of C(ti,t3) are shown 
in Fig. Q, where the variance anisotropy ratio for datasets B and D is plotted as a function of X4(C(ii,i3)) for the 
values C(ti,t^) ~ 0.23, 0.33 and 0.44. This figure shows that the variance anisotropy ratio grows with the strength 
of the dynamical heterogeneity, for all three fixed values of (7(^1,^3). 

According to the hypothesis, dynamical heterogeneity originates in the Goldstone modes associated to fluctuations 
in the time reparametrization (see Eq. Q). Therefore, the hypothesis implies that the correlation length of the 
dynamical heterogeneity should be similar to the correlation lengths of the transverse components of the variables tti 
and 7r2, and the longitudinal variable a should be less correlated in space. Our results, shown in Fig. (|8]), show that 
this is indeed the case. The correlation volume corresponding to the transverse fluctuations, V"^ , closely tracks the 
one corresponding to the dynamical heterogeneities, Vc^^ , and they both grow together as the temperature is reduced 
or the timescale is increased. By contrast, for the longitudinal fluctuations a we find that their correlation volume is 
small and essentially constant; it does not correlate with the correlation volume of the dynamical heterogeneity, nor 
with the temperature or the time scale. 
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FIG. 4: 25% probability contours for dataset B with C{ti,t3) = 0.23. Contour B corresponds to the 25% of the 
probabihty density shown in Fig. Contour B' corresponds to the same times as contour B, but with a much 
smaller coarse graining size. Contour label B" corresponds to the same coarse graining size as B, but with much 

shorter times. 




FIG. 5: Left panel: Ratio between the variances of transverse and longitudinal fluctuations as a function of the 
initial time ti. Right panel: Ratio between the correlation volumes of transverse and longitudinal fluctuations as a 
function of ti. The color coding identifies the datasets and is the same as in Fig. (IT|). 
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FIG. 6: Left panel: Ratio between the variances of transverse and longitudinal fluctuations as a function of the 
strength of the fluctuations of (7(^1,^3), measured by X4 [C(<i, ^3)]. Right panel: Ratio between the correlation 
volumes of transverse and longitudinal fluctuations as a function of Xi [C'l^ii^a)]- The color coding is the same as in 

Fig. (g. 
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FIG. 7: Ratio between the variances of transverse and longitudinal fluctuations as a function of X4 [^{^1,^3)], for 
systems B and D for the values (7(^1,^3) = 0.23, 0.33, and 0.44. The color coding is the same as in Fig. (fll). 
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FIG. 8; Correlation volumes for transverse and longitudinal fluctuations, plotted against the correlation volume for 

Ci3 = (7(^1,^3). The color coding is the same as in Fig. (fll). 



IV. CONCLUSIONS 

In this paper we tested the hypothesis that dynamical heterogeneity arises from Goldstone modes related to a 
broken continuous symmetry under time reparametrizations. In other words, we tested whether dynamical hetero- 
geneity is associated with the presence of spatially correlated fluctuations in the time variables. To verify this, we 
studied probability distributions that allowed us to distinguish between time reparametrization fluctuations (trans- 
verse fluctuations) and other fluctuations (longitudinal fluctuations) . We also tested for possible correlations of both 
the strength and the correlation volume of the fluctuations in the time variable with the dynamical susceptibility X4, 
which is normally used to probe dynamical heterogeneity. Altogether, we found that at the lowest temperatures, for 
the longest timescales and for the largest coarse graining lengths, the transverse fluctuations became stronger than 
the longitudinal fluctuations, which is consistent with the hypothesis, and that this was the case independently of 
the system considered. We also found that the correlation volumes of the time reparametrization fluctuations were 
proportional to the correlation volume of the dynamical heterogeneity, while the correlation volumes of the longitu- 
dinal fluctuations were small and independent of the correlation volumes of the dynamical heterogeneity. All of this 
was despite the fact that, to simplify the analysis, we imposed some extra conditions on the form of the correlations, 
which may have made the agreement with the hypothesis appear less good than it would have been otherwise. 

With regards to universality, the evidence we found was mixed. On the one hand, there were clear differences 
in the details of the results for different systems, for example for the anisotropy ratios. On the other hand, all the 
trends we observed were the same across systems, and the results for the correlation volumes did show some hints 
of universality, although the relatively large noise in this measurement did not allow for definite conclusions to be 
drawn. In any case, the commonality in the results is strong enough to suggest that other systems may display similar 
qualitative behaviors. Thus, we expect that it would be very instructive to apply the same kind of test to data from 
other slowly relaxing systems, such as particle tracking data from glassy colloids [5^12j and from granular systems 
close to jamming [l4 | [T5 ] . 

Finally, considering the success of the tests presented here, it becomes natural to ask if it is possible to extract from 
the data the actual fluctuating reparametrization 4>p{t), and to study its properties directly. In fact, Ref. (36j already 
shows some progress in that direction. 
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